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We study the charge screening effect in the hadron-quark mixed phase. By including the charge 
screening effect, rearrangement of charged particles occurs and some part becomes locally charge- 
neutral. As a result the equation of state for the mixed phase becomes close to that given by the 
Maxwell construction, which means that the Maxwell construction would effectively gain the 
physical meaning again even in the system with two or more chemical potentials. We also discuss 
the interplay of the surface tension and the Coulomb interaction. Both effects would restrict the 
region of the mixed phase in the core of hybrid stars. 
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1. Introduction 

It is widely believed that quark matter exists in high-temperature and/or high-density situa- 
tions like in the relativistic heavy-ion collisions (RHIC) §j§ or in the core of neutron stars [Q, g|, 
and the "deconfinement transition" has been actively searched. Theoretical studies using model 
calculations or based on the first principle, lattice QCD have been also carried out by many 
authors to find the critical temperature. Although many exciting results have been reported, the 
deconfinement transition has not been clearly understood yet. 

We, hereafter, consider the phase transition from hadron phase to three-flavor quark phase in 
neutron-star matter in a semi-phenomenological way. Since many theoretical calculations have 
suggested that the deconfinement transition should be of first order in low-temperature and high- 
density area ^], we assume the first-order phase transition in this paper. If the deconfinement 
transition is of first order, one may expect a mixed phase during the transition. Actually, the hadron- 
quark mixed phase has been considered during the hadronization era in RHIC [J7|, ||, |9p or in the 



core region of neutron stars [|10|, |llj, |12fl . 

There is an issue about the mixed phase for the first-order phase transitions with two or more 
chemical potentials Jl3|]. We often use the Maxwell construction (MC) to derive the equation 
of state (EOS) in thermodynamic equilibrium, as in the liquid-vapor phase transition of water, 
where both phases consist of single particle species. However, if many particle species participate 
in the phase transition as in neutron-star matter, MC is no more an appropriate method. Before 



Glendenning first pointed out [13], many people have applied MC to get EOS for the first-order 
phase transitions [14, 15, [T^, 17] expected in neutron stars, such as pion or kaon condensation and 
the deconfinement transition. 

Let us consider a deconfinement transition in neutron-star matter. We must introduce many 
chemical potentials for particle species, but the independent ones in this case are reduced to two, i.e. 
baryon-number chemical potential /Ib and charge chemical potential jUq, due to chemical equilib- 
rium and total charge neutrality. They are nothing but the neutron and electron chemical potentials, 
\l n and n e , respectively. In the mixed-phase these chemical potentials should be spatially constant. 
When we naively apply MC to get EOS in thermodynamic equilibrium, we immediately notice 
that /Ib is constant in the mixed-phase, while jU e is different in each phase because of the difference 
of the electron number density in these phases. This is because MC uses EOS of bulk matter in 
each phase, which is of locally charge-neutral and uniform matter; many electrons are needed in 
hadron matter to compensate the positive charge of protons, while in quark matter charge neutrality 
is almost fulfilled without electrons. Thus 

in MC, where superscripts "Q" and "H" denote the quark and hadron phases, respectively. Glen- 
denning emphasized that we must use the Gibbs conditions (GC) in this case instead of MC, which 
relaxes the charge-neutrality condition to be globally satisfied as a whole, not locally in each phase 



P13|] . GC imposes the following conditions, 

M B Q = 

pQ = pH T Q = T H (1 2 ) 
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He demonstrated a wide region of the mixed phase, where two phases have a net charge but totally 
charge-neutral: EOS thus obtained, different from that given by MC, never exhibits a constant- 
pressure region. He simply considered a mixed phase consisting of two bulk matters separated by 
a sharp boundary without any surface tension and the Coulomb interaction, which we call "bulk 
Gibbs" for convenience. "Bulk Gibbs" requires that each matter can have a net charge but the total 
charge is neutral, 



0, 



(1.3) 



/vP c Q h + (Wv)p c H h 

O H 

where f v means the volume fraction of quark matter in the mixed phase and "pY' " means charge 
density in each matter. Figure [l] shows the phase diagram in the jUb - He plane. We can see that 
there is a discontinuous jump in \i e for the case of MC, while the curve given by "bulk Gibbs" 
smoothly connects uniform hadron matter and uniform quark matter; the mixed phase can appear 



in a wide jUb region in "bulk Gibbs", in contrast with MC [18J. 
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Figure 1: Phase diagram in the — fi e plane. There appears no region of the mixed phase by the calculation 
with MC, while a wide region of the mixed phase by the "bulk Gibbs" calculation. 



However, this "bulk Gibbs" is too simple to study the mixed phase, since we must consider 
matter with non-uniform structures instead of two bulk uniform matters; the mixed phase should 
have various geometrical structures where both the number and charge densities are no more uni- 
form. Then we have to take into account finite-size effects like the surface and the Coulomb 
interaction energies. We show that the mixed phase should be narrow in the /Ib space by the charge 
screening effect, and derive EOS for the deconfinement transition in neutron-star matter. We shall 
see EOS results in being similar to that given by MC. We also discuss the interplay of the Coulomb 
interaction effect and the surface effect in the context of the hadron-quark mixed phase. Preliminary 
results for the droplet case has been already reported in Ref. [|9j|. 

2. Brief review of the previous works 

Heiselberg et al. p0| ] studied a geometrical structure in the mixed phase by considering the 
spherical quark droplets embedded in hadron matter. They introduced the surface tension and 
treated its strength as a free parameter because the surface tension at the hadron-quark interface 
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has not been clearly understood. They pointed out that if the surface tension parameter a is large 
(a > 90 MeV/fm 2 ), the region of the mixed phase is largely limited or cannot exist. Subsequently 
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Figure 2: Schematic image of structured mixed phase. 



Glendenning and Pei [10] have suggested "crystalline structures of the mixed phase" which have 
some geometrical structures, "droplet", "rod", "slab", "tube", and "bubble" (Fig. ||), assuming a 
small q[|l0j |2lj ]. The finite-size effects are obvious in these calculations by observing energies. We 
may consider only a single cell, by dividing the whole space into equivalent Wigner-Seitz cells: the 
cell size is denoted by Ryy and the size of the lump (droplet, rod, slab, tube or bubble) by R (Fig. 
[3]). Then the surface energy density is expressed in terms of the surface tension parameter a as 



R ' 



(2.1) 



where d denotes the dimensionality of each geometrical structure; d = 3 for droplet and bubble, 
d = 2 for rod and tube, and d = 1 for slab. The Coulomb energy density reads 



£c = 27re 2 (p c H h -p c Q h ) R 2 ® d (f v 



<Wv) 



2(d-2)- [ 1 



+ fv 



(d + 2) 



-1 



(2.2) 
(2.3) 



where we simply assume an uniform density distribution in each phase. When we minimize the 
sum of £5 and Sc with respect to the size R for a given volume fraction fy, we can get the well- 
known relation, 

£ 5 = 2£ C . (2.4) 

This implies that an optimal size of the lump is determined by the balance of these finite-size 
effects. Eventually we can express the sum of the surface and Coulomb energy densities with an 
optimal cell size 



3f v d 



2d 



*d(fv) 



1/3 



(2.5) 
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Thus we can calculate the energy of any geometrically suuctured mixed phase with (2.5) by chang- 
ing the parameter d. Many authors have taken this treatment for the mixed phase[ 10, 18[ 2(J. Note 
that the energy sum in Eq. ( jl^l ) becomes larger as the surface tension gets stronger, while the 



nucleon phase 
quark phase 



relation Eq.(|2.4|) is always kept. 

However, this treatment is not a self-consistent but a perturba- 
tive one, since the charge screening effect for the Coulomb potential 
or the rearrangement effect of charged-particles in the presence of the 
Coulomb interaction is completely discarded. We shall see that the 
Coulomb potential is never weak in the mixed phase, and thereby this 
treatment overestimates the Coulomb energy. The charge screening ef- 
fect is included only if we introduce the Coulomb potential and con- 
sistently solve the Poisson equation with other equations of motion for 
charged particles. Consequently it is a highly non-perturbative effect. 
Norsen and Reddy [ p2| ] have studied the Debye screening effect in the 
context of kaon condensation to see a large change of the charged- 
particle densities like kaons and protons. Maruyama et al. have numer- 
ically studied it in the context of liquid-gas phase transition at subnu- 




clear densities [|23[], where nuclear pasta can be regarded as geometrical 
structures in the mixed phase. Subsequently, we have also studied kaon 



nucleon quark 



Figure 3: Wigner-Seitz 
approximation for the 
droplet case. R is the 
droplet radius and /?w 
the cell radius. 



condensation at high-densities [24], where kaonic pasta structures ap- 
peared in the mixed phase. Through these works we have figured out 
the role of the Debye screening in the mixed phase. 



Voskresensky et al. [25] explicitly studied the Debye screening ef- 
fect for a few geometrical structures of the hadron-quark mixed phase. 
They have shown that the optimal value of the size of the structure can- 
not be obtained due to the charge screening, if the surface tension is 
not so small. They called it mechanical instability. It occurs because 
the Coulomb energy density is suppressed for the size larger than the Debye screening length (cf. 



Eq. ( |4.2[ )). They also suggested that the properties of the mixed phase become very similar to those 
given by MC, if the charge screening effectively works. The apparent violation of the Gibbs con- 
ditions in MC (Eq. flLll )) can be remedied by including the Coulomb potential in a gauge-invariant 
way: the number density of charged particles is given by a gauge-invariant combination of the 
chemical potential and the Coulomb potential, and thereby it can be different in each phase for 
a constant charge chemical potential if the the Coulomb potential takes different values in both 
phases. However, they used a linear approximation to solve the Poisson equation analytically. 

If the Coulomb interaction effect is so large, it would be important to study it without recourse 
to any approximation. In this paper we numerically study the charge screening effect on the struc- 
tured mixed phase during the deconfinement transition in neutron-star matter in a self-consistent 
way. Actually we shall see importance of non-linear effects included in the Poisson equation. 



3. Self-consistent calculation with the Coulomb potential 



We have presented our formulation in detail in Ref. [26, 27]. Here we only briefly explain it. 
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We consider the geometrically structured mixed phase (SMP) where one phase is embedded in the 
other phase with a certain geometrical shape. We divide the whole space into equivalent charge- 
neutral Wigner-Seitz cells. The Wigner-Seitz cell is a representative of equally divided space under 
approximation of a geometrical symmetry: sphere in three dimensional (3D) case, rod in 2D case 
and slab in ID case with the size R^. The cell consists of two phases in equilibrium: one phase 
with a size R is embedded in the other phase as illustrated in Fig. [3| 

Quark phase consists of u, d, s quarks and electron. We incorporate the MIT Bag model and 
assume the sharp boundary at the hadron-quark interface: u and d quarks are treated to be massless 
and s is massive (m s = 150MeV), and they interact with each other by the one-gluon-exchange 
interaction inside the bag. Hadron phase consists of proton, neutron and electron. The effective 
potential is used to describe the interaction between nucleons and to reproduce the saturation prop- 
erties of nuclear matter. 

Then total thermodynamic potential (Qtotai) consists of the hadron, quark and electron contri- 
butions and the surface contribution: 

n to tal=n H + nQ + £ls, (3.1) 

where ^h(q) denotes the contribution of hadron(quark) phase. We here introduce the surface con- 
tribution Q.s, parametrized by the surface tension parameter a, Q.s = oS with 5 being the area of the 
interface. Note that it may be closely related with the confining mechanism and unfortunately we 
have no definite idea about how to incorporate it. Actually many authors have treated its strength 
as a free parameter and seen how the results are changed by its value[|l(j, 18, EO], which we also 
follow in this report. 

To write down the thermodynamic potential we use the idea of the density functional theory 
(DFT) within the local density approximation (LDA) [ 28[ 29|] . The Coulomb interaction energy is 



included in £2h(q) an d can be expressed in terms of particle densities, 



E V = 1 -Y( d\ d y Q ^) Q ?f\ (3.2) 
2fjJv w \r-r'\ 



where i = u,d,s,p,n,e with £?, being the particle charge (Q = — e < for the electron) and V w the 
volume of Wigner-Seitz cell. Accordingly the Coulomb potential is defined as 

V(r) = -Z[ dV e -pn + V , (3.3) 

where Vb is an arbitrary constant representing the gauge degree of freedom. We fix the gauge by a 
condition V(Rw) = in this paper. Operating a Laplacian V 2 on the Coulomb potential V(r), we 
automatically derive the Poisson equation. 

There are introduced six chemical potentials in Eq. ( |3.1[ ) corresponding to particle species 
in the cell. Then we have to determine now eight variables, i.e., six chemical potentials, jUi (i = 
u,d,s,p,n,e), and the radii R and Ryy. Fixing R and Rw, we have four conditions due to chemical 
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equilibrium in each phase and at the interface: 

Md = Ms, 

fa = fa + 2fa, 

(^ = 2^ + ^). (3.4) 

Therefore, once two chemical potentials fa and fa are given, we can determine other four chemical 
potentials, fa, fa, fa and fa. We determine fa by the global charge neutrality condition: 

/vP c Q h + (Wv)p c H h = 0, (3.5) 

where the volume fraction fv = (j^J > an d d denotes the dimensionality of each geometrical 
structure. At this point fy is still fixed. The particle densities p,- are given by the equations of 
motion, which are derived form the extremum condition, = 0, 

^^kin+str ,~ ^, 

where £kin+str stands for the kinetic and strong interaction energy except the Coulomb interaction 
energy. 

Then we find the optimal value of 7? (Rw is fixed and thereby fy is changed by 7?) by using 
one of GC; 

pQ =p n + Pa (37) 

where the pressure coming from the surface tension is given by 

dS 



dV Q 



(3.8) 



The pressure in each phase pQ( H ) j s given by the thermodynamic relation: pQ( H ) = — £2q( H )/Vq( H ), 
where £2q(h) is the thermodynamic potential in each phase and given by adding electron and the 



Coulomb interaction contributions to £2q(h) in Eqs. ( |3.1| ). Finally, we determine R w by minimizing 
thermodynamic potential. Thus, once /Ib is given, all other values fa (i = u,d,s,p,e) and R, R w 
can be obtained. Accordingly, the particle density profiles are given by Eq. ([3i]). 

We numerically solve Eqs. ( |3.6| ) and the Poisson equation in a self-consistent way. In nu- 
merical calculation, every point inside the cell is represented by a grid point (the number of grid 
points /Vgnd ~ 100). Equations of motion are solved by a relaxation method for a given baryon- 
number chemical potential under constraints of the global charge neutrality. Note that we keep GC 
throughout the numerical procedure. 

4. Numerical results and discussions 

First, we present the thermodynamic potential density of each uniform matter and that given 
by "bulk Gibbs" in Fig. ^ In uniform matter, hadron phase is thermodynamically favorable for 



7 



Charge screening effect on hadron-quark mixed phase in compact stars 



Tomoki Endo 




Mb < 1225 MeV and quark phase above it. Therefore we plot 8(0, difference of the thermodynamic 
potential density between the mixed phase and each uniform matter: 

f ffltotai - ^ niform Mb < 1225MeV, 
8(o={ ^ H (4.1) 

1 Ototai - «Q niform Mb > 1225MeV, 

where (Ot oU i = n tota i/W, etc. 

The results are shown in Figs. |5| and |6[ In these figures we also depict two results for com- 
parison: one is given by"bulk Gibbs", where the finite-size effects are completely discarded. The 
other is the thermodynamic potential given by a perturbative treatment of the Coulomb interac- 
tion, which is denoted by "no screening"; discarding the Coulomb potential V(r), we solve the 
equations of motion to get uniform density profiles (no rearrangement), then evaluate the Coulomb 
interaction energy by using the density profiles thus determined. Note that "no screening" violates 



the gauge invariance, while our treatment respects it []25|]. Since chemical potential itself is gauge 
variant, we have to include V(r) together like in Eq. to satisfy the gauge invariance: when the 
Coulomb potential is shifted by a constant value, Vif) ==> Vir) — Vq, the charge chemical potential 
should be also shifted as Mi Mr +NjVo. Incidentally, the phase diagram in the Mb — Me plane 
(see, e.g., Fig. |l|) is not well-defined by this reason, while many people have written it [18, ^(J. We 



can see the screening effects by comparing the results given by "no screening" with those by the 
self-consistent calculation denoted by "screening". 8(0 given by MC appears as a point denoted 
by a circle in Figs. || and ^ where two conditions, P Q = P H and Mb = Mb > WQ satisfied. On the 
other hand the mixed phase derived from "bulk Gibbs" appears in a wide region of Mb- Therefore, 
if the region of the mixed phase becomes narrower, it signals that the properties of the mixed phase 
become close to those given by MC. One may clearly see that ft)t ota i becomes close to that given by 
MC due to the finite-size effects, the effects of the surface tension and the Coulomb interaction. 

The large increase of 8(0 from the "bulk Gibbs" curve comes from the effects of the surface 
tension and the Coulomb potential. Since the surface tension parameter is introduced by hand, 
we must carefully study the effects of the surface tension and the Coulomb interaction, separately. 
From the difference between the result given by "no screening" and that by "bulk Gibbs", we can 
roughly say that about 2/3 of the increase comes from the effect of the surface tension and 1/3 
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Figure 5: Difference of the thermodynamic po- 
tential density as a function of baryon-number 
chemical potential fi B for rj = 40 Me V/fm 2 . If 8 CO 
is negative, the mixed phase is a thermodynami- 
cally favorable state. MC determines one point of 
phase transition in uniform matter, denoted as a 
circle in the Hb-8(0 plane. 
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Figure 6: Same as Fig. | for C7 = 60 Me V/fm 2 . 
The negative 8(0 region is narrower than the a = 
40 Me V/fm 2 case. 



from the Coulomb interaction (see Eq. (2A)). Comparing the result of self-consistent calculation 
with that of "no screening", we can see that the change of energy caused by the screening effect is 
not so large , but still the same order of magnitude as that given by the surface effect. 
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Figure 7: Density profiles in the droplet phase 
given by "no screening" for = 1 189 MeV and 
a = 60 Me V/fm 2 . They are uniform in each 
phase. R = 7.2 fm and Rw = 12.8 fm. 
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Figure 8: Density profiles and the Coulomb po- 
tential given by the self-consistent calculation for 
the same parameter set as Fig. 5. R = 7.7 fm and 
R w = 18.9 fm. 



If the surface tension is stronger, the relative importance of the charge screening effect be- 
comes smaller and the effect of the surface tension becomes more prominent, as is seen in Figs. || 
andg. 

Although the charge screening has not so large effects on bulk properties of the matter, we 
shall see that it is remarkable for the charged particles to change the properties of the mixed phase. 
The screening effect induces the rearrangement of the charged particles. We can see this effect 
by comparing Fig. [7| with Fig. [|. The quark phase is negatively charged and the hadron phase 
is positively charged. The negatively charged particles in the quark phase such as d, s, e and 
the positively charged particle in the hadron phase p are attracted toward the boundary. On the 
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contrary the positively charged particle in the quark phase u and negatively charged particle in the 
hadron phase e are repelled from the boundary The charge screening effect also reduces the net 
charge in each phase. In Fig. ||, we show the local charge densities of the two cases shown in 
Figs. f7] and || The change of the number of charged particles due to the screening is as follows: 
In the quark phase, the numbers of d and s quarks and electrons decrease, while the number of u 
quark increases. In the hadron phase, on the other hand, the proton number should decrease and 
the electron number should increase. Consequently the local charge decreases in the both phases. 
In Fig. ^ we can see that the core region of 
the droplet tends to be charge-neutral and near 
the boundary of the Wigner-Seitz cell is almost 
charge-neutral. In Figs. [K] and 1 1 we present the 
lump and cell radii for each density. As we have 



shown in the previous paper Q25|], the Coulomb 
energy is suppressed for larger R by the screen- 
ing effect. The R dependence of the total ther- 
modynamic potential comes from the contribu- 
tions of the surface tension and the Coulomb 
interaction: the optimal radius giving the mini- 
mum of the thermodynamic potential is then de- 
termined by the balance between two contribu- 
tions, since the former gives a decreasing func- 
tion, while the latter an increasing one. If the 
Coulomb energy is suppressed, the minimum of 
the thermodynamic potential is shifted to larger 
radius. As a result the size of the embedded 
phase (R) and the cell size (Rw) become large. 
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Figure 9: Local-charge densities for the "no 
screening" case and the case of the self-consistent 
with the screening effect. For "no screening", 
charge density of each phase is constant over the 
region. The absolute value of the charge density 
is larger than that given by the self-consistent cal- 
culation in each phase. In the hadron phase, the 
charge density becomes almost vanished near the 
cell boundary r = R\y. 

In Ref . [|25j] they demonstrated that the minimum 

disappears for a large value of the surface tension parameter: the structure becomes mechanically 
unstable in this case. 

We cannot show it directly in our framework because such unstable solutions are automatically 



excluded during the numerical procedure, while we can see its tendency in Figs. 10 and 11: R and 
Rw get larger by the screening effect. 

We also see the relation between the size of the geometrical structure and the Debye screening 
length. The Debye screening length appears in the linearized Poisson equation and is then given as 



-2 



An^Qf 

f 



'd(pf) 



dfif 



-2 



AkQ p 



dfip 



(4.2) 



where (pf 1 ) stands for the averaged density in quark phase, (Pp h ) is proton number averaged den- 
sity in the hadron phase and (p t ? h ) is the electron charge density averaged inside the cell. It gives 
a rough measure for the screening effect: At a distance larger than the Debye screening length, the 
Coulomb interaction is effectively suppressed. 



In Fig. 10 we show sizes of geometrical structure for "no screening" case. If we ignore the 
screening effect, the size of the embedded phase is comparable or smaller than the corresponding 
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Figure 10: Lump and cell radii given by the 
"no screening" calculation. The Debye screening 
length is also depicted for comparison. R is thick- 
solid line and Rw is thick-dashed line. We can see 
the size of the structure becomes less than the De- 
bye screening length. 
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Figure 11: Same as Fig. |l0| given by the self- 
consistent calculation with the screening effect. 
The size of the structure becomes larger than that 
given by "no screening", and consequently ex- 
ceeds the Debye screening length. 




Figure 12: Schematic graphs of the droplet size and the Debye screening length. Right figure shows the 
case of the self-consistent calculation with the screening effect and left figure "no screening". 



quark Debye screening length Ap (Fig. 12). This may mean that the Debye screening is not so 
important. Actually, many authors have neglected the screening effect due to this argument^, ^(]] . 
In Fig. [ll], however, we see that the size of the embedded phase can be larger than (Fig. 12) in 
the self-consistent calculation. We can also see the similar situation about Rw and Xfy. This means 
that the screening has important effects in this mixed phase. We cannot expect such a effect without 
solving the Poisson equation because of the non-linearity. 

We present EOS in Figs. 13 and 14. We can see that EOS of two cases become close to that 
given by MC by including the finite size effects. Moreover, EOS of "screening" becomes more 
similar to that given by MC than that of "no screening". We show the pressure results in Figs. [15] 
and 16, which is given by P = Pb^§^~ ov P = —D./V: the former is given by using Figs. |13| and 
14 and the latter by using the results of Fig. |5[ The pressure becomes more similar to that given by 
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MC by the charge screening effect, which shows a larger pressure near the beginning and weaker 
one near the termination of the phase transition. 
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Figure 13: Equation of state for cr = 40 MeV/fm 2 
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Figure 16: Same as Fig. |15| given by the self- 
consistent calculation with the screening effect. 



We have used the surface tension parameter a in the present study. Surface tension is a very 
difficult problem because it should be self-consistent with the two phases of matter, quark and 
hadron. Lattice QCD, based on the first principle, predicts that o may be 10-100 MeV/fm 2 [|3~l|, 
Although this result is for high temperature case, our choice is within it. Moreover, the results given 
by other model calculations [|3|, 34 , |35 , are similar to our value. 

Let us briefly consider some implication of these our results for compact star phenomena. 



Glendenning QlQp suggested many SMP appear in the core region by using "bulk Gibbs": the mixed 
phase should appear for several kilometers. However we can say that the region of SMP should be 
narrow in the jUb space and EOS is more similar to that of MC due to the finite-size effects. These 



results look to be consistent with those given by other studies. Bejger et al. [37] have examined the 
relation between the mixed phase and glitch phenomena, and shown that the mixed phase should 
be narrow if the glitch is generated by the mixed phase in the inner core. On the other hand the 



gravitational wave asks for density discontinuity in the core region[38]. 
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5. Summary and concluding remarks 



We have studied the charge screening effect in the hadron-quark mixed phase in neutron-star 
matter. We have elucidated the charge screening effect comparing the results of "no screening" 
with those of "screening" which includes non-linearity of the Poisson equation. The density pro- 
files of the charged particles are drastically modified by the screening effect, while the thermody- 
namic potential is not affected so much; the charge rearrangement induced by the screening effect 
tends to make the net charge smaller in each phase. Consequently the system tends to be locally 
charge-neutral, which suggests that the Maxwell construction (MC) is effectively justified even if 
it is thermodynamically incorrect. In this context, it would be interesting to refer to the work by 
Heiselberg [p9|], who studied the screening effect on a quark droplet (strangelet) in the vacuum, 
and suggested the importance of the rearrangement of charged particles. 

We have seen that thermodynamic quantities such as thermodynamic potential and energy 
become close to those derived from MC by the screening effect, which also suggests that MC is 
effectively justified due to the screening effect. As another case of the system with two or more 
chemical potentials, kaon condensation has been also studied [23] and the results are shown to 
be similar to those in the present study. Thus the importance of the screening effect should be a 
common feature for the first-order phase transitions in high-density matter. 

We have included the surface tension at the hadron-quark interface, while its definite value is 
not clear at present. An uncertainty in the value of the surface tension does not allow to conclude 
whether the mixed phase exists or not. There are many estimations for the surface tension at 
the hadron-quark interface in lattice QCD [ [31] , 32], in shell-model calculations [33, [34], 35] and 
in model calculations based on the dual-Ginzburg Landau theory [p6|]. Our parameter is in that 
reasonable range. 

We have considered some implications of our results for neutron star phenomena. The screen- 
ing effect would restrict the allowed region of the structured mixed phase in neutron stars, in con- 
trast with a wide region given by "bulk Gibbs"[10, 21]. It could be said that they should change 
the bulk property of neutron stars, especially the structure of the core region. 

Compact stars have the strong magnetic field and its origin is not well understood. One possi- 
bility is that it comes from ferromagnetism of the quark matter in the core ©Ullll. Therefore it 
should be interesting to include the magnetic field contribution in our formalism. We have assumed 
zero temperature here. It would be much interesting to include the finite-temperature effect. Then 
it is possible to draw the phase diagram in the He - T plane and we can study the properties of 
the deconfinement phase transition; our study may be extended to treat the mixed phase to appear 
during the hadronization of QGP in the nucleus-nucleus collisions and supernova explosions. 

In this study we have used a simple model for quark matter to figure out the finite-size effects 
in SMP. However, it has been suggested that the color superconductivity would be a ground state 



of quark matter []T8|, [43], B4fl. Hence we will include it in a further study. The hadron phase should 
be also treated more realistically; for example, we should include the hyperons or kaons in hadron 
matter. In the recent studies the mixed phase has been also studied |4^, 47] in the context of 
various phases in the color superconducting phase. 
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